Anomalous criticality in the quantum Hall transition at n = Landau level of graphene 

with chiral-symmetric disorders 
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We investigate numerically whether the chiral symmetry is the sole factor dominating the crit- 
icality of the quantum Hall transitions in disordered graphene. When the disorder respects the 
chiral symmetry, the plateau-to-plateau transition at the n — Landau level is shown to become 
anomalous (i.e., step-function-like). Surprisingly, however, the anomaly is robust against the in- 
clusion of the uniform next-nearest neighbor hopping, which degrades the chiral symmetry of the 
lattice models. We have also shown that the ac (optical) Hall conductivity exhibits a robust plateau 
structure when the disorder respects the chiral symmetry. 

PACS numbers: 73.43.-f, 72.10.-d, 71.23.-k 
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I. INTRODUCTION 



Kicked off by the recent observation of the characteris- 
tic quantum Hall effect in graphene^— the effect of disor- 
der has attracted considerable attention4rZ A key inter- 
est that is specific to the electronic structure of graphene 
is the chiral (A-B sub-lattice) symmetry^ and the scat- 
tering between two Dirac cones that reside at K and K' 
points in k-space. Accordingly, the nature of disorder 
(i.e., whether the disorder preserves the chiral symmetry 
or not) should be essential in discussing the electronic 
properties of disordered graphene. A potential disorder 
breaks the chiral symmetry, whereas a disorder in bonds 
preserves the symmetry. The spatial range of disorder is 
another important factor, since the inter- valley scattering 
occurs when the disorder is short-ranged while strongly 
suppressed when long-ranged. 

In our previous papers^, we have considered the disor- 
der in the hopping energies as a model for a ripple ; 10 i 11 
which is an intrinsic disorder in graphene. We have 
clearly demonstrated that for such a chiral-symmetry 
preserving disorder, the criticality of the quantum Hall 
transition at the n = Landau level becomes anoma- 
lously sensitive to the spatial correlation of disorder: 
When the spatial range of the bond disorder exceeds 
a few lattice constant, no broadening of the graphene- 
specific n — Landau level occurs and the associated 
step-function like Hall transition exhibits an anomalous 
(almost fixed-point) criticality. This sharply contrasts 
with the case of potential disorder which degrades the 
chiral symmetry, where the criticality is similar to the 
ordinary quantum Hall transition^ 

Now, when one describes the actual graphene with 
a tight-binding model, the effect of the next-nearest 
neighbor hopping in the honeycomb lattice is not 
negligibleJ^r— Since the next-nearest neighbor hopping 
t' connects sites in the same sub-lattice, the chiral sym- 
metry of the original the honeycomb lattice is broken, so 
that one might naively think that the speciality about 



the n = Landau level associated with the chiral sym- 
metry should be washed away. It is to be noted, how- 
ever, that the space group (i.e., rotation and reflection 
symmetries as well as the translation symmetry) of the 
honeycomb lattice are preserved, and, indeed, the two 
Dirac cones still exist at K and K' points at an energy 
shifted from zeroJ£ The existence of the doubled Dirac 
cones has a topological origin. For instance, they exist 
even in the presence of a third- neighbor hopping t" in one 
of the three directions (which is equivalent to the it flux 
model) with t' = as shown in Ref. H, where the lattice 
continuously crosses over to the square lattice as t" is 
increased. In such a system, the Dirac cones are located 
at general points away from K and K' in the Brillouin 
zone. If we add if in this case, tilted Dirac cones appear 
at a non-zero energy. The tilted Dirac cones are interest- 
ing, since they have been observed in an organic system 
a-(BEDT-TTF) 2 I 3 under a high pressureJI^ 

Given these variety of situations, we want to clarify 
in the present paper whether the anomalous criticality 
is preserved or washed out. First, we shall clarify the 
question: Do we have to discriminate the presence or 
otherwise of the chiral symmetry in the underlying hon- 
eycomb lattice (before we turn on the disorder) and the 
presence or otherwise of the chiral symmetry in the dis- 
order? In order to study these, we adopt the honey- 
comb lattice model rather than the effective Dirac model, 
since the chiral symmetry has to do with the AB sublat- 
tices. To examine the robustness of the criticality for the 
n = Landau level, we calculate the Hall conductivity in 
the honeycomb lattice as well as the it flux model with 
the bond disorder to examine the influence of the next- 
nearest neighbor hopping. We find that the anomalous 
criticality at the n = Landau level for the long-ranged 
bond disorder is robust against the inclusion of the uni- 
form next-nearest neighbor hopping, namely the break 
down of the chiral symmetry of the lattice models, which 
yields, in general, the energy shift as well as the tilt of the 
Dirac cones. The Hall conductivity for non-interacting 
electron systems is mathematically equivalent to a topo- 
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logical integer called a Chern number^ - — Taking into 
account the recent progress in numerical evaluation of the 
Chern number,— ~— we adopt a formulation based on the 
lattice gauge technique to evaluate the Hall conductivity 
in disordered finite systems. The formulation is based on 
the topological invariant used in the lattice gauge the- 
ory, which can also be understood as a two dimensional 
generalization of the King-Smith- Vanderbilt formula for 
polarizationi 24 i 26 The present formulation has been suc- 
cessfully applied to the calculation of the Hall conductiv- 
ity not only for uniform systems but also for disordered 
systems M- 

Our second question is: How the criticality appears in 
the ac (optical) Hall conductivity, a xy (uj)? The optical 
response of QHE system is now attracted growing inter- 
est. For instance, the cyclotron emission from graphene 
QHE system has been studied theoretically^ as well as 
experimentally with infrared transmission. 29 Specifically, 
Morimoto et &)M have shown that the optical Hall con- 
ductivity <J xy {ijj) has a Hall plateau structure even in the 
ac (~ THz) regime, both in the ordinary two-dimensional 
electron gas (2DEG) and in graphene in the quantum 
Hall regime, although the plateau height is no longer 
quantized in ac. In graphene (Jxyiyi), which reflects the 
unusual Landau level structure, the plateau structure re- 
mains against a significant strength of disorder due to an 
effect of localization. For 2DEG the THz spectroscopy 
has recently been performed, and a plateau structure of 
optical Hall conductivity has been detected in the THz 
region. 31 Hence it is interesting to investigate the con- 
sequence of the chiral symmetry in the ac Hall effect in 
graphene. 



II. CHIRAL SYMMETRY IN LATTICE MODELS 

Let us begin with the chiral symmetry. The symme- 
try is defined by the existence of a local operator 7 with 
7 2 = 1 that anti-commutes with the Hamiltonian H as 
{7, H} = 0. For bipartite lattices such as honeycomb and 
square, the lattice can be decomposed into A and B sub- 
lattices. If we take 7 = riieA ex P( 17rc l Ci )' ^ ne f erirnon 
operator ci is transformed as 7^7 = — Cj for i G A and 
jdj = a for i G B. We can readily show that the lattice 
model that has only hopping terms iycjcj connecting A 
and B sites has the chiral symmetry {7, H} = 0, irre- 
spective of the details of the hopping terms ty . 

An important consequence of the chiral symmetry is 
that the energy levels appear in pairs, {E,—E}, and 
hence the energy spectrum is symmetric around E = 0. 
Assume that ipE is an eigenstate with energy E satisfying 
HipE = EipE- Then the state jipE is an eigenstate with 
energy —E, since H{^ipE) = —jHipE = —E(j^e)- For 
E = the eigenstates tpE=o and 7?Ab=o ar e degenerated. 
We can rearrange the eigenstates into ijj± = (1 ± 7)?/>£ = o 
to make them simultaneous eigenstates of 7 with "/ip± = 
±V>±- This implies that the system has an additional 
symmetry in the zero-energy space. From the definition 



of 7, it is clear that the eigenstate VM-(-) nas its ampli- 
tudes only on the B(A) sub-lattice sites. The criticality 
of the Hall transition at E = can thus be sensitive to 
the presence or absence of the chiral symmetry^ 



III. DC HALL CONDUCTIVITY 

A. Hamiltonian 

The tight-binding Hamiltonian for the honeycomb lat- 
tice with the next-nearest neighbor hopping is given by 

nn nnn 

H = J2 kje-^clcj + £ t'er^dc, + h.c, 

where the summation is performed for nearest neighbor 
sites in the first term and for next-nearest neighbor sites 
in the second term (Fig. [T]) . Here iy and t' are assumed 
to be real, while the Peierls phases, {%} and are 
determined such that the sum of the phases around a 
loop is equal to the magnetic flux piercing the loop in 
units of the flux quantum O = h/e. We specify the 
strength of the uniform magnetic field by the magnetic 
flux <fi per hexagon in the honeycomb lattice. We intro- 
duce a bond randomness in the nearest-neighbor transfer 
energy as ty = — £+£iy, where t and SUj denote the uni- 
form and the disordered components, respectively. The 
disorder in 5tij is assumed to have a Gaussian distribu- 
tion P{5t) = e- <5 * 2 / 2o ' 2 /V27rcr 2 with a variance a. The 
spatial correlation length 77 in the random components is 
specified by requiring 

{5U 3 5t kl ) = {5f)e-^- r ^/^\ 

where ry represents the bond iy, and () the ensemble 
average^ 3 -. The length is hereafter measured in units of 
the distance a between nearest-neighbor sites. It is to 
be noted that the Hamiltonian H is chiral symmetric 
if t' = even in the presence of disorder in the nearest 
neighbor hopping ty , while the chiral symmetry is broken 
for t' 0. In actual graphene, the magnitude of t' is 
estimated to be 0.2 - 0.02tH 



B. Hall conductivity as a Chern number 

In order to evaluate the Hall conductivity, let us 
consider a L x x L y system, and impose boundary 
conditions for the single-particle wave function \ip) as 
T x ( y )(L x(y) )\if) = e^v^ip), where T x{y) (L x{y) ) denotes 
the single-particle magnetic translation operator in the 
x(y) direction. The string gauge in the corresponding 
brick-layer lattice is adopted to treat this twisted bound- 
ary conditions. 34 When the energy gap exists above the 
Fermi energy Ee, the Hall conductivity related 
to a topological integer Ce called the Chern number 
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FIG. 1: Left: Honeycomb lattice with nearest neighbor hop- 
ping —t (solid lines) and next-nearest neighbor hopping t' 
(dotted lines). Right: Energy dispersion for <f> = 0,t' /t — 
0.15. Dirac cones are located at K — (2ir/3-\/3a, 2-7r/3a) and 
K' = (47r/3^3a,0) with energy -3t' M 
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where denotes the wave function of the Fermi sea in 
the twisted boundary conditions specified by the phase 
= (</'xj0y);(O < < j ) x(y) < 27r) with the normalization 
(0|0) = 1 and T 2 is the two-dimensional parameter space 
for the phase 0. The Chern number can be expressed, by 
introducing a S(j) x by 5(j) y mesh for the parameter space 
T 2 , as C E = ^ F m<n with 



F m , n = <b A-d<t>, A = (ip\V$ip). 

JdS mn 

Here S mn stands for a 8(j) x by 8<f> y square of the mesh 
with points = (mS(j) x ,n54>y), ((m + l)S(j> x , nS<f> y ), ((m + 
l)6<f> x , (n + l)8(j)y), (mS(j) x , (n + l)S<p y ). The surface in- 
tegral on S mn is replaced by the line integral by Stokes' 
theorem. It is useful to note that this line integral F m ^ n 
can be written as 

F m , n = In 14(0) [7,(0 + A<f> x )U- 1 (<f>+ A4> y )U;\4>), 

with U a {4>) = exp(fjp + ^ < ^ a A a d(f> a ) (a — x,y), where 

A4> x = {6<f> x ,Q), A(f) y = (0,6</> y ) and 0= (m54> x ,n5<j) y ). 

Following Fukui et al.j2i we consider a link variable 
defined on the mesh as 

- (0(0)|0(0 + A0J) 

1(0(0)10(0 + A0J)| 

and define the Chern number Ce on the mesh as Ce = 
SiE ro /m,n, where 

F m ,„ = In U x (0) 0, (0 + A0 X ) L>- 1 (0 + A0. y ) £7- 1 (0) 

with — 7T < F m>n /\ < 7r. The Chern number Cg defined in 
such a way on the mesh reduces to the Chern number Ce 
defined on the continuum space in the limit 5<f> x ( y ) — > 0. 
It is to be noted that Ce is gauge invariant, and takes 



an integer value for an arbitrary mesh width 5(j> x r y y We 
have then Ce = Ce for a fine enough meshi 24 ' 25 

In the present case, the wave function of the Fermi 
sea |0) is a Slater determinant of the single-particle 
wave function \ip), so that the link variable U a (<fi) 
for an M particle state is given by 2 ^ U a (<fi) = 
det F a (0)/| det Vq(0)|, where V a {<j>) is an M x M ma- 
trix. The elements of V are defined by (V a (4>))ij = 
(^(0)|^(0 + A0J), (1 < i,j < M)j, with the i-th 
single-particle eigenstate \ipi(4>)) satisfying H\ipi((j>)) = 
£i\Wi{(f>)) (ei < ... < £fc < £fc+i,...) in the boundary 
condition specified by 0, where the Fermi energy Ee is 
assumed to lie between Sk and £k+i- The energy gap at 
the Fermi energy is mostly guaranteed for finite and ran- 
dom systems due to the level repulsion. We thus adopt 
this formulation to evaluate the Chern numbers for each 
sample. By taking an average over samples in each en- 
ergy bin, we obtain ensemble-averaged Hall conductivity 
as a function of the Fermi energy. 

C. Anomalous criticality in the n — level 

The Chern number Ce for the case of a bond- 
disordered honeycomb lattice with t' = 0, where the chi- 
ral symmetry is exactly preserved, is shown in Figj^J We 
can clearly see that the plateau transition in the n = 
Landau level at E = becomes almost step-function like 
as a function of the Fermi energy as soon as the spatial 
correlation length, 77, of the bond disorder exceeds the 
lattice constant a, as previously shown in Ref. |7|. Re- 
markably, this behavior is fixed-point like in that it shows 
almost no system-size dependence. The anomalous Hall 
transition is accompanied by the absence of broadening 
of the n = Landau level. 

Next we show what happens when we switch on the 
next-nearest neighbor hopping t' 7^ 0. Note that the 
Dirac cones still exist at K and K' at an energy, — 3t', 
shifted from zeroJ^ The Hall conductivity in the case of 
t' /t = 0.05 is shown in Fig. [3]for several values of the dis- 
order correlation length 77. It is clearly seen that the Hall 
plateau transition at the n = Landau level is again 
anomalously sensitive to the spatial correlation of the 
bond disorder. The results for different system sizes are 
displayed in Fig(U where we sec little size-dependence at 
n = 0. This occurs despite the fact that the chiral sym- 
metry of the lattice model is broken by the next-nearest 
neighbor hopping t' ^ 0. Consequently, the energy spec- 
trum is not symmetric about E = with the n = 
Landau level shifted from E = 0. If we vary the strength 
of the next-nearest neighbor hopping t 1 , the critical en- 
ergy for the anomalous Hall transition associated with 
the n — Landau level is simply shifted with t' (FigO. 
For n = ±1 Landau levels, by contrast, the usual scaling 
behavior (i.e., a narrower transition region for a larger 
system) is observed (Fig. 2]). 

The present results amount to that the anomalous 
transition remains at the energy of the Dirac cones, where 
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FIG. 2: (Color online) The Chern number Ce (Hall conduc- 
tivity in units of e 2 /h) as a function of the Fermi energy E/t 
for various values of the spatial correlation length rj/a of the 
bond disorder, for t'/t — 0, (f>/4>o = 1/50, a/t = 0.115 and 
L x /(VSa/2) = L v /(3a/2) = 20. 



FIG. 4: (Color online) The Chern number Ce (Hall con- 
ductivity in units of e 2 /h) as a function of the Fermi en- 
ergy E/t in the presence of the next-nearest-neighbor hopping 
t'/t = 0.05 for rj/a = 1.5, for cj>/(j>Q = 1/36, a/t = 0.115 and 
(L x /(VSa/2), L y /(Sa/2)) = (24, 18) and (12, 12). 
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FIG. 3: (Color online) The Chern number Ce (Hall conduc- 
tivity in units of e 2 /h) as a function of the Fermi energy E/t in 
the presence of the next-nearest-neighbor hopping t'/t = 0.05 
for rj/a = 0.0,0.5, 1.0 and 1.5, for </>/</> = 1/36, a/t = 0.115 
and (L x /(V3a/2),L y /(3a/2)) = (24, 18). 



FIG. 5: (Color online) The Chern nmber Ce (Hall conduc- 
tivity in units of e 2 /h) as a function of the Fermi energy 
E/t for various values of the next-nearest-neighbor hopping 
t', with a fixed rj/a = 1.5, <j>/(j>o = 1/36, a/t — 0.115 and 
{L x /(V3a/2),L y /(3a/2)) = (24,18). 



the breakdown of the usual chiral symmetry does not af- 
fect the anomalous criticality although the energy spec- 
trum is affected. This stability of the anomalous critical- 
ity is attributed to the preserved effective chiral symme- 
try for the 2x2 low-energy effective Hamiltonian in the 
vicinity of the Dirac cones 



D. Tilted Dirac cones 

To examine further the robustness of the anomaly at 
the n = Landau level, we consider a lattice model hav- 
ing an third-neighbor hopping t" across the hexagon of 



the honeycomb lattice (Fig. [H]). It has been shown® that, 
for the case of t' = 0, this model has also two Dirac 
cones in the first Brillouin zone for 1 < i" < 3, which 
are adiabatically connected to the Dirac cone in the orig- 
inal honeycomb lattice as the hopping f is continuously 
varied. When t" is increased to t (with t' = 0), the model 
is equivalent to the n flux model, which has been ana- 
lyzed as a typical model with the Dirac dispersion.— If 
we further introduce the second-neighbor hopping t' ^ 
in the model, we have not only shift the energy of the 
Dirac cones, but also tilt the cones. For simplicity, here 
we confine ourselves to the case of t" — t, where two 
anisotropic Dirac cones are located at (7r/2\/3a, 7t/2a) 
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FIG. 6: The n flux model having additional hopping t ' (thick 
dashed lines) in the y direction across the hexagon of the 
honeycomb lattice. 
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FIG. 7: (Color online) Tilted Dirac cones for the 7r-flux model 
with t" = t and t'/t = 0.15. The Dirac cones are located at 
(k x ,k y ) = (-K/2y/3a,Tv/2a) and (— iv/2y/Sa, 5iv/6a) and E — 
-0.3t. 



and (— 7r/2\/3a, 57r/6a). In the presence of t' , Dirac cones 
are located at E = — It' . 

The Hall conductivity obtained for this model is shown 
in Fig. [8] in the presence of spatially correlated disorder 
with r]/a = 1.5. For the case oit' = 0, the system is chiral 
symmetric and two Dirac cones exist at E = 0. We thus 
obtain anomalous criticality at E = 0. When the hopping 
t' is introduced, thereby degrading the chiral symmetry, 
we clearly see that the anomalous plateau transition are 
shifted in energy but preserved. The shift of the critical 
energy with t' is again consistent with the energy shift of 
the Dirac cones. The result suggests that the anomaly in 
the n = Landau level is preserved even when the Dirac 
cones are tilted, although the tilt angle achieved in the 
present calculation is rather small (up to ~ 5°). 
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FIG. 8: (Color online) The Chern number Ce (Hall conduc- 
tivity in units of e 2 /h) as a function of the Fermi energy E/t 
in the presence of the next nearest-neighbor hopping t' for 
the 7r-flux model, for r\/a = 1.5, 4>/4>o = 1/36, cr/t = 0.115 
and {L x /(VSa/2),L y /(3a/2)) = (18, 12). The strength of the 
hopping t" across the hexagon is assumed to be the average 
of \tij\ around the corresponding hexagon. 



IV. AC HALL EFFECT 

Now we would like to look at the ac extension of the 
Hall conductivity, i.e. optical Hall conductivity <j xy (ui) 
to discuss the effect of the chiral symmetry. The crucial 
interests is whether the existence or otherwise of the chi- 
ral symmetry dominates ac responses as well, where we 
pay particular attention to the plateau structure in the 
optical Hall conductivity a xy (ui). 

So we have calculated a xy (uj), where we solved the 
Hamiltonian containing various disorders with the exact 
diagonalization method, since we are interested in the 
localization effects in the ac Hall conductivity. We then 
compute the optical Hall conductivity by using the Kubo 
formula, since the method for calculating the topological 
(integer) Chern numbers in the dc case is not applicable 
to the optical conductivity. For the ac Hall conductivity, 
we consider the usual honeycomb lattice with different 
types of disorder to examine how the symmetry of disor- 
der is reflected in the ac conductivity. 

The Kubo formula for the ac conductivity is given as 




(1) 



where e a is the eigenenergy, j x the current matrix 
elements between the eigenstates, and Ef the Fermi 
energy2£. The current operator is given by a gauge field 
derivative of the Hamiltonian as 

. _ dH 
3 ~ dA' 
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FIG. 9: (a) a X y(eF, ui) plotted against Fermi energy ef and 
frequency uj for the honeycomb tight-binding model with a 
potential disorder yj {Se 2 } = O.lt which breaks the chiral 
symmetry, (b) g xv (ef,uj) for uj — 0,0.5t/h is displayed to 
show the broadening of the step structure associated with 
the n — Landau level. The parameters are 4>/4>o = 1/18, 
(L x /(y/Sa/2), L v /(3o/2)) = (27, 24). 

where the gauge field and the Peierls phase in the Hamil- 
tonians are related via 2ir8ij — A ■ (r, — rj). 

To see the effect of the chiral symmetry on a xy (oj), we 
have examined two cases: (i) The system with a site- 
potential disorder introduced as 

H = ]T -te- 2 ^c\ Cj + h.c. + fc * c K 

with a random potential Se having a Gaussian distri- 
bution. The potential disorder, which breaks the chi- 
ral symmetry, models an effect of charged impurities in 
graphene samples, (ii) The system with random hopping, 

H = ^2(-t+ %)e -2,rie «cjcj + h.c, 

which preserves the chiral symmetry, is considered to be 
coming from ripples in graphene samples. So we would 
look into the effect of these on the optical Hall conduc- 
tivity. 



Fig[9]shows <T xy (eF, w) plotted against Fermi energy ef 
and frequency u> with a potential disorder yj (Se 2 ) = O.lt 
which breaks the chiral symmetry. The curve for to = 
in FiglHJa) corresponds to dc Hall conductivity. If we 
vary Fermi energy ef for a fixed frequency ui, we no- 
tice that there are two distinct regions. For the Dirac- 
cone-like energy region bounded by the van Hove singu- 
larities (\ef\ < t) we observe step structures at Dirac 
Landau levels, where we have wider plateaus for higher 
Landau levels due to a y/B dependence of cyclotron en- 
ergy for Dirac QHE. At the van Hove singularities the 
Dirac QHE crosses over to the ordinary QHE for usual 
fermions around the band edges, and outside the van 
Hove singularities (\ef\ > t) an usual QHE step sequence 
is recognized. 

If we look at the dependence on the frequency cj, each 
plateau exhibits the cyclotron resonance behavior as uj 
is increased, where the resonance shape differs between 
the two regions. For the usual QHE region the resonance 
occurs for small ui because of small cyclotron frequency 
(oc B). In the Dirac QHE region, we have resonances for 
larger u (oc y/B), which appear at various resonance fre- 
quencies that correspond to resonances between different 
for LL's with non-uniform LL spacings with a peculiar 
selection rule, \n\ O \n + 1| for Dirac QHE. An inter- 
esting point observed for the honeycomb tight-binding 
model is that these two regions emerge from the nature 
of the honeycomb band dispersion and separated by van 
Hove singularities. 

FigEJb) depicts the step structure around n = LL 
for dc and ac (uj — O.ht/h) Hall conductivities, which 
shows that the n = steps are blurred both for dc and 
ac response when the disorder does not respect the chiral 
symmetry. 

Now we come to the original question of what happens 
to the ac Hall conductivity when the disorder respects the 
chiral symmetry. Figure flOl shows a xy (EF,uj) for the case 
of a random hopping yj {St 2 ) — O.lt which respects the 
chiral symmetry. The overall structure of the (J X y(^F, 
is similar to the random potential case, but we can imme- 
diately notice that the n = step in a xy (EF,uj) exhibits 
an anomalously sharp step structure at even for the op- 
tical Hall conductivity as in the dc Hall conductivity. 
Indeed, FigfTOTa) indicates that a xy (sF — 0, ui) shows a 
sharp step for all the values of w, which should come from 
the protected n — LL by the chiral symmetry. 

We can conclude that the chiral symmetry and the zero 
mode protection give a crucial effect on the ac response 
of the system as in the dc response, where the presence 
of the chiral symmetry gives a rise to the anomalously 
sharp step structure in <j xv (ef,u>) associated with delta- 
function like n = LL. The honeycomb lattice calcu- 
lation, performed here for the optical Hall conductivity, 
is not only more realistic than the effective Dirac field 
description, but also reveals an existence of two regions 
(Dirac QHE and usual fermionic QHE) coming from the 
van Hove singularity in the band dispersion, where totally 
different plateau width and cyclotron resonance struc- 
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FIG. 10: (a) a xy (eF,u) plotted against Fermi energy ef 
and frequency cj for the honeycomb tight-binding model with 
random hopping v/ {St 2 ) = O.lt which respects the chiral sym- 
metry, (b) a xy (EF,u) with uj = 0,0.5t/fi is displayed to 
show the anomalous step structure at n=0 LL retained for 
ac Hall conductivity. The parameters involved are r\ja — 2.0, 
<j>/<f>o = 1/18, (L x /(VSa/2),L y /(3a/2)) = (27,24). 



anomalous criticality associated with the 5-function like 
density of states at the n = Landau level in graphene, 
which yields a step-function like Hall plateau transition. 
The ac as well as dc Hall conductivity has been exam- 
ined as a function of the Fermi energy in the presence 
of the spatially-correlated chiral symmetric bond disor- 
der. For the dc Hall conductivity, it has been shown that 
the anomalous criticality is surprisingly robust against 
the introduction of the next-nearest neighbor hopping in 
the honeycomb lattice, which is considered to be appre- 
ciable in real graphene. The anomalous criticality may 
therefore be realized in a suspended clean graphene where 
ripples are the major disorder. As for the symmetry, the 
present results implies that the anomalous criticality can 
be insensitive to the breakdown of the usual chiral sym- 
metry as long as the effective chiral symmetry that pro- 
duces the massless Dirac cones is preserved. We have also 
shown that the anomaly is robust against a slight tilt in 
the Dirac cones, which suggests the possibility to observe 
the anomalous criticality not only in graphene but also in 
organic materials. For the ac Hall conductivity, we have 
shown that even for finite frequencies, the transition at 
the n = Landau level exhibits the anomalous critical- 
ity, with the robust plateau structure retained for long- 
ranged, chiral-symmetric bond disorders. Based on these 
numerical results, we conclude that the chiral symmetry 
manifests itself both in the dc and the ac Hall transitions 
yielding the anomalous criticality at the n = Landau 
level of graphene. 
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